Spectral properties of a partially spin-polarized one-dimensional Hubbard/Luttinger superfluid 
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We calculate the excitation spectra of a spin-polarized Hubbard chain away from half-filling, using a high- 
precision momentum-resolved time-dependent Density Matrix Renormalization Group method. Focusing on the 
U <Q case, we present in some detail the single-fermion, pair, density and spin spectra, and discuss how spin- 
charge separation is altered for this system. The pair spectra show a quasi-condensate at a nonzero momentum 
proportional to the polarization, as expected for this Fulde-Ferrel-Larkin-Ovchinnikov-like superfluid. 

PACS numbers: 



Systems of interacting fermions obeying Fermi liquid the- 
ory exhibit a one-to-one correspondence between their low- 
energy quasiparticle excitations and those of a non-interacting 
Fermi gas. The quasiparticles have renormalized energy and 
spectral weight, but possess the same charge and spin quan- 
tum numbers as the corresponding noninteracting fermions. 
This scenario breaks down in one dimension (ID): in this case, 
each Fermi surface reduces to two points in momentum space, 
at fc = ±:kp. The resulting Fermi-surface nesting is present at 
all densities and spin polarizations and destabilizes the Fermi 
liquid, converting it instead to a Luttinger liquid ||il|2j,|3|], even 
for a weak interaction. 

In a Luttinger liquid with zero spin polarization, the el- 
ementary excitations are collective density fluctuations that 
carry only either spin ("spinons") or charge ("holons"). These 
excitations have different dispersions, and, obviously, do 
not carry the same quantum numbers as the original "bare" 
fermions. This leads to the spin-charge separation picture, in 
which a fermion injected into the system separates ("fraction- 
ates") into an (anti-)holon and a spinon, each of them carrying 
a share of the fermion's quantum numbers. The phenomenon 
of spin-charge separation, and fractionation of particles more 
generally, is an important concept in strongly-correlated sys- 
tems, and has intrigued physicists for decades. Its signatures 
have been observed experimentally in ID organic conductors 
(j4|], metallic wires carbon nanotubes and nanowires 
in semiconducting heterostructures [FtI]. Proposals have been 
made to seek for evidence of these phenomena in cold atomic 
gases iSil. 

In the ID Hubbard model, the low-energy spin and charge 
modes of the Luttinger liquid decouple as long as the system 
either is at half-filling or has zero spin-polarization. However, 
if the system is away from half-filling and has a nonzero mag- 
netization, the collective modes that constitute the elementary 
low-lying excitations are linear combinations of the spin and 
charge fields fT3, U , 12], so although one still has a Luttinger 
liquid with fractionalized fermions, it is no longer strictly a 
"spin-charge separation" scenario. The field-theoretical for- 
mulation of the Luttinger liquid theory has proven very ef- 
fective in describing the low-energy physics of a variety of 
models. However, a fully quantitative and general picture of 
how the spin and charge degrees of freedom couple to form 



full-fledged fermions is still missing. 

In this work, we will focus our attention on the negative- 
U (attractive) ID Hubbard model, away from half-filling and 
at nonzero spin polarization. This model can now be studied 
experimentally with ultracold atoms in an optical lattice. As 
usual, the Hamiltonian is 
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where c^^ creates a fermion with spin a = f, J, at site £; 
nia — cl^Cia', t is the hopping matrix element, which we set 
to unity (we also set the lattice spacing to unity); and U is the 
interaction strength that in this work will be considered nega- 
tive (attractive). The negative and positive U versions of this 
model can be mapped exactly onto each other by the "canon- 
ical" transformation that applies a particle-hole and momen- 
tum change to one spin species. Thus our results are general, 
and can be translated to the positive-t/ case [13]. 

For large negative U the fermions form tightly bound pairs 
that behave as hard-core bosons flS*] . These bosons are pre- 
vented from fully condensing in ID due to quantum fluctua- 
tions. They form a "quasi-condensate", with pair correlations 
that decay as a power-law that, in some regime of parame- 
ters (large \U\ > At) can dominate over the single-fermion 
correlations at large distances lll4ll . In the polarized case, the 
ground-state of this system is the ID version of the Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) superfluid Q [H [H, 
in which the pairs forming the quasi-condensate have nonzero 
center of mass momentum ±Q with Q = kp^ — kpi, where 
kpa is the Fermi momentum of the fermions with spin a. This 
was confirmed numerically in Ref. jlTll . and subsequent stud- 
ies Hill). 

The Luttinger-liquid and FFLO aspects of this system can 
be heuristically understood as follows: At large negative U, 
the spin-polarized ground state consists of empty sites (O's), 
sites occupied by pairs (2's) and excess up fermions (fs), with 
sites singly-occupied by J,'s being only "virtual" states. The 
density of excess f s is Q/tt. An "f exchanges positions with 
the O's and 2's with hopping t and thus moves with bandwidth 
4t. At half-filling, the background the | moves through is half 
O's and half 2's, so the relative motion moves spin but no den- 



sity on average: this collective mode is then purely a spinon. 
But if we move away from half-filling, the number of O's dif- 
fers from the number of 2's, so when an t moves, it on average 
moves some density as well as spin: this light (bandwidth At) 
mode of the Luttinger liquid is then not purely spin, but in- 
stead is a particular linear combination of spin and charge (we 
will call this light mode "spinon-like")- In the limits of nearly 
complete polarization or either zero or complete filling, the 
t's become just regular fermions carrying the full charge and 
spin . This scenario has been confirmed numerically in Ref. 
1 12, 20i 21(1 . by looking at the real-time evolution of spin and 
charge distributions. 

At large negative U the 2's do not move freely past the 
O's; this exchange happens via a virtual intermediate un- 
paired state with energy \U\, resulting in effective hopping 
teff = —2t^/U. Thus this motion of 2's relative to O's con- 
stitutes the heavy "holon-like" mode of the Luttinger liquid 
with a smaller bandwidth. Also, when a 2 moves past an t, 
the ground state has a sign change. This means the wave- 
function of the quasi-condensate of bosonic 2's has a node 
at each If these nodes were equally-spaced, this would 
be an FFLO standing-wave condensate with momentum ±Q. 
However, the f 's actually form a ID Luttinger liquid with di- 
vergent position fluctuations, so the momentum distribution 
of the pairs instead has a power-law divergence at ±Q\ this 
ID partially spin-polarized superfluid state should perhaps be 
termed "quasi-FFLO". 

The Hamiltonian ([T]l can be solved exactly by means of the 
Bethe Ansatz 122. .23., .24.1 , and the dispersion of the elemen- 
tary excitations can be obtained |25, 3. 27]. However, the 



20 

1.6 



actual Green's functions and spectral properties can only be 
calculated in certain limits [28,1, and numerical methods have 
been crucial to fill in the blanks and compare to experiments 
[n the following, we use the time-dependent exten- 
sion of the Density Matrix Renormalization Group (tDMRG) 
E. 32] method to obtain estimates for various Green func- 
tions in real-time and real-space with unprecedented accuracy 
1 331] ■ To extract the dynamical response of the system, we cal- 
culate the correlators G{x ~x' ,t' -t) = i{0{x',t')0^ {x, <)), 
where O is an operator of interest. Fourier transforming then 
yields the correspondin g sp ectral weights as functions of mo- 
mentum and frequency ll3lLl33Ll3411 : 
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where Eq is the ground state energy, and the sum runs over all 
the eigenstates of the system, with energy £"„ . All the results 
will be plotted using a log-scale for the intensity, with several 
orders of magnitude between the intensities of the weakest 
and strongest features. At very small scales, some ripples or 
oscillations appear as a consequence of the numerical Fourier 
transform, and the commensuration of the lattice. These ef- 
fects get amplified near zero momentum and frequency. 

In Fig. [TJa) and (b) we show the dynamic structure factor 
for the charge and spin densities, respectively, for an unpolar- 
ized Hubbard chain at quarter-filling (in this paper we always 
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FIG. 1: (color online) Dynamical structure factors of the (a) 
"charge" density n{k, oj), and (b) spin Sz{k, to) for an unpolarized, 
quarter-filled Hubbard chain with U = —St. (c) Spectral weights for 
adding {lo > /j.) or removing (u < fermion for the same system; 
/i is the chemical potential. Frequencies are in units of the hopping 
t = 1. The colors are set by the logarithm of the spectral intensity. 



use L — 80 and U — —8t). The charge excitations display 
gapless modes at momenta fc = and k = ±2kF = ±7r/2, 
and a continuum ranging from uj^Otouj = t = 4ie//. 
This spectrum is formed primarily by holon-antiholon excita- 
tions. It is qualitatively similar to the particle-hole spectrum of 
the corresponding non-interacting system, but with a reduced 
bandwidth. However, this system is a superfluid with a spin 
gap of = 5t, as is seen in the spectral weight of the spin (Fig. 
lb); this is the energy "cost" of breaking a Cooper pair. The 
spinon has band-width = At, and the spectral weight of Sz 
vanishes strongly as fc — > 0, since the total spin is conserved 
and the matrix element for making spin excitations thus van- 
ishes at zero momentum. 

The single-particle spectral weight for the quarter-filled, 
unpolarized system is shown in Fig. [Uc), where we plot the 
imaginary part of the one-particle Green's function. The up- 
per and lower features, for positive and negative frequencies, 
correspond to the inverse photoemission (IPES) and photoe- 
mission (PES) spectra, resulting from adding or removing a 
fermion, respectively. We have shifted the energies relative 
to flie chemical potential /i = [Eo{N + 1) - Eo{N - l)]/2, 
which lies in the center of the gap for this unpolarized system. 
This gap is a manifestation of the spin gap due to Cooper pair- 
ing: the ground state is a total spin singlet with all fermions 
paired. The added fermion has no "partner" to pair with, while 
removing a fermion requires breaking an existing pair, so both 
processes are gapped. 

Again, we can heuristically understand many features of 
these spectra using the large-negative-J7 description discussed 
above. The unpolarized ground state is a quasi-condensate 
of 2's that form a Luttinger liquid of repulsively interacting 
bosons. An added f forms a spinon and much of its spectral 



3 




-0,8 -0.4 0,0 0.4 0.8 -0.8 -0.4 0.0 0.4 0.8 



FIG. 2: (color online) "Photoemission" spectra for quarter-filled 
spin-polarized Hubbard chain, with U = —8t, = 24, A*'; — 16. 
The energy scale has been chosen relative to the average chemical 
potential /i — {jj,^ + /i|)/2. The Fermi levels for each spin species 
are indicated by the horizontal lines (see text). 



weight thus follows a spinon dispersion with bandwidth At. 
Since the wavefunction changes sign when the | exchanges 
position with a 2, the lowest-energy spinon states are at the 
momenta ±7r/4 set by the density of the 2's. However, the 
added fermion may also excite holon modes, and a careful 
look at the upper part of Fig. [TJc) reveals a continuum, with 
a weaker feature at the lower edge of the continuum which 
has a holon-like dispersion. This continuum arises when part 
of the added momentum is used to excite holon modes of the 
quasi-condensate. 

Removing a fermion requires breaking a pair (a "2"), and 
this process apparently couples more strongly to the holon de- 
grees of freedom, as can be seen by the flatter dispersion of 
the strong part of the spectrum at low momentum in the bot- 
tom part of Fig. [Tfc). However, at higher momentum, this 
PES spectrum, although much weaker, has a continuum with 
a mostly spinon-like dispersion. Here the process apparently 
removes a low-momentum pair from the quasi-condensate and 
makes a spinon, with the spinon taking most of the momen- 
tum. At half-filling, there is particle-hole symmetry and the 
PES and IPES spectra are thus equivalent, both containing 
strong spinon and holon signals. ll29ll 

We now turn our attention to the single-fermion spectrum in 
the polarized case, shown in Fig. [21 where we took N^—N[ ~ 
8. Since the system is no longer symmetric under time- 
reversal, the spectral functions for the up and down fermions 
are different. Correspondingly, we can determine the chemi- 
cal potential of each species: ji^r = [E{Na- + 1) — E{Nc — 
l)]/2. The average chemical potential is /i = + 
while the effective Zeeman field is /i = {jj.-^ — iJii)/2. We 
find it instructive to plot the spectra with energies relative to 
ji. Note that /i is still in the "pairing" gap, but now /i| is in the 
band above the gap, while is in the band below the gap. 



At this fairly large \U\, we can describe this system as a 
quasi-condensate of bosonic 2's with density kpi/T^ = 1/5 
and momentum ±Q = ivr/lO, plus a density Q/tt = 1/10 of 
excess unpaired f's. The 2's are bound pairs and sit below the 
gap and just below /ij^. In the PES spectrum one can remove 
a fermion of either spin from one of these singlet pairs; these 
are the strong low-momentum features near uj — ji = —3. 
The weaker bands dispersing strongly to lower energy from 
these features arise from removing one member of a pair and 
leaving the other member in a spinon-like state. 

The excess unpaired |'s lie at energy just below /i| and can 
be seen there in the spin-up PES spectrum. The wavefunctions 
of the t's change sign on passing each 2; as a result the lowest 
energy states of the corresponding spinon-like modes are at 
±fci?|; it is near these momenta where the up PES intensity 
is largest. The strongest bands in the up spectrum cross 
at -izkp^, just as in the noninteracting system. But one can 
also see weaker bands crossing at ±(2fci?| — kp]), which 
coiTespond to three-particle excitations in the noninteracting 
system. At low momentum and energies below /i| there is 
a fairly flat holon-like dispersion of the up spectral weight; 
presumably here the excitation also transfers some momentum 
to the (heavy) pairs. 

In the spin-down IPES spectrum, there is a heavy holon- 
like band at energies just above This arises from adding 
a down fermion that pairs with one of the excess up fermions 
with momentum |fc| < kp^, resulting in a pair (a 2) which car- 
ries most of the added momentum. This feature in the IPES 
is strong only for |fc| > kpi, since the down-spin states at 
lower momentum than this are already occupied. This band 
continues to the PES spectrum below crossing the chemi- 
cal potential at A: = ±fcFX, as in the noninteracting system. At 
the zone boundary this holon-like band splits in to two faint 
features at {oj — fi) near -1 and -1.5, for reasons we do not 
yet understand. In the PES spectrum there are also weaker 
features approaching /i| at momenta ±(2A;f| — kpi), which 
again correspond to three-particle excitations in the noninter- 
acting system. 

In the IPES spectrum above the gap, we can see that 
for both spins the added particle can excite a continuum 
of states with both spinon-like and holon-like dispersions. 
The sharpest feature is a spinon band, which is substantially 
sharper in the up IPES than in the down. The general appear- 
ance of this part of the IPES spectral weight is similar to that 
of the unpolarized case in Fig. 1(c). 

In order to examine the quasi-condensate and the pair-like 
excitations, we calculated the spectral weight of removing a 
pair hi — Ci|Ci|, shown in Fig. [Sjc) and (d). For the unpolar- 
ized case, we see a large concentration of spectral weight at 
zero energy and zero momentum in the PES, due to the zero- 
momentum quasi-condensate. In the polarized case, the zero- 
energy spectral weight of this quasi-condensate splits into two 
features at fc = ±Q = ±(fc_F| — kpi), as expected for this 
FFLO-like state. 

In Fig. Oa) we show the dynamical structure factor for the 
density operator in the polarized case. One noteworthy differ- 
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FIG. 3: (color online) Dynamical structure factors of the (a) 
"charge" density n{k,Lj), and (b) spin lowering operator S-{k,ij) 
for the polarized, quarter-filled Hubbard chain of Fig. 2. (c) Spectral 
weight for removing a pair from the unpolarized system, and (d) the 
polarized system. In (c) and (d) the energy scale is relative to the 
chemical potential of a pair 2fi. 



ence from the unpolarized case is the appearance of a weak 
spinon-like feature. This occurs because in this polarized sys- 
tem away from half-filling, the spinon-like mode is no longer 
purely spin so couples to the density. 

Since time-reversal symmetry is broken by the spin polar- 
ization, the response functions for the spin operators Sz, S+, 
and S- are now all different. The structure factor for Sz (not 
shown) exhibits, besides the excitations across the gap present 
in the unpolarized system (Fig. 1(b)), a gapless band that 
originates from spin excitations within the bands that cross 
the Fermi surfaces. The operator 5*+ flips spins up, breaking 
pairs, so its spectrum only shows excitations across the gap. 
The richest of these spin dynamical functions is S-{k,uj), 
shown in FiglSjb). The action of the operator S- on the 
ground state can cause three possible outcomes that each oc- 
cupy a separate energy window: 1) it can break a pair; 2) it 
can flip an unpaired t to an unpaired |, making a gapless spin 
fluctuation; or 3) the flipped spin can pair with another un- 
paired t and be absorbed by the quasi-condensate. We find 
that the spectral weight for this last process is very weak, but 
detectable; it is visible just below energy -4 in Fig. 3b. We be- 
lieve this very small weight is due to the product of two small 
factors: the low probability (due to fermionic antisymmetry) 
that two unpaired fs are on adjacent lattice sites before one 
of them is flipped down, and the low overlap between the re- 
sulting state after flipping, and the ground state with a bound 
pair, since the latter mostly consists of doubly-occupied sites. 

To summarize, we have reported and discussed the rich fea- 
tures of the particle, pair, spin and density spectral weights for 
the quasi-FFLO superfluid ground state of a partially spin- 
polarized fermionic Hubbard chain with attractive interac- 



tions. We have found that a rigorous treatment, particularly 
to describe properties involving excitations, should still rely 
on the Luttinger liquid picture. 
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